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We propose a scheme for subwavelength localization of an atom conditioned upon the absorption 
of a weak probe field at a particular frequency. Manipulating atom-field interaction on a certain 
transition by applying drive fields on nearby coupled transitions leads to interesting effects in the 
absorption spectrum of the weak probe field. We exploit this fact and employ a four-level system 
with three driving fields and a weak probe field, where one of the drive fields is a standing-wave 
field of a cavity. We show that the position of an atom along this standing wave is determined when 
probe field absorption is measured. We find that absorption of the weak probe field at a certain 
frequency leads to subwavelength localization of the atom in either of the two half- wavelength regions 
of the cavity field by appropriate choice of the system parameters. We term this result as sub-half- 
wavelength localization to contrast it with the usual atom localization result of four peaks spread 
over one wavelength of the standing wave. We observe two localization peaks in either of the two 
half-wavelength regions along the cavity axis. 



I. INTRODUCTION 

Precision measurement of the position of an atom pass- 
ing through a standing-wave field has attracted consid- 
erable attention in recent years. Several schemes have 
been proposed for the localization of an atom using opti- 
cal methods 0. High-resolution position measurements 
of the atom with optical techniques are of considerable 
interest from both a theoretical and an experimental 
point of view. Interest in the area is largely due to its 
applications to many areas of optical manipulations of 
atomic degrees of freedom, such as laser cooling QjBose- 
Einstein condensation 0, and atom lithography Q and 
measurement of the center-of-mass wavefunction of mov- 
ing atoms 

It is well known that optical methods provide better 
spatial resolution in position measurements of atoms. For 
example, in the optical virtual slits scheme the atom in- 
teracts with a standing-wave field and imparts a phase 
shift to the field. Measurement of this phase shift then 
gives the position information of the atom [||. Another 
related idea based on phase quadrature measurement is 
considered in Ref. Q ■ Kunze et al. Q demonstrated how 
the entanglement between the atomic position and its 
internal state allows one to localize the atom without di- 
rectly affecting its spatial wave function. It is shown that, 
by using Ramsey interferometry, the use of a coherent- 
state cavity field is better than the classical field to get a 
higher resolution in position information of the atom Q . 
Resonance imaging methods have been used in experi- 
mental studies of the precision position measurement of 
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the moving atoms Q, 0, 0] . 

More recently, atom-localization methods based on the 
detection of the spontaneously emitted photon during the 
interaction of an atom with the classical standing-wave 
field are considered 0, 0, 0, 0] . We consider some 
of these proposals in detail to contrast them with the 
current proposal. Qamar et al. 14] suggested a simple 
scheme for localization of an atom by using a simple two- 
level system interacting with the classical standing-wave 
field. They showed that the frequency of the sponta- 
neously emitted photon carries information about the po- 
sition of the atom. In another related study, they used a 
three-level atom, where the upper two levels are driven by 
a classical standing-wave field and spontaneously emitted 
photon measurement from the upper level to lower level 
gives information about the atomic position |l5j . It has 
been shown that coherent control of spontaneous emis- 
sion in multilevel system gives line-narrowing and even 
spontaneous emission quenching |l6j . By using three 
drive fields, the phase and amplitude control of the driv- 
ing field on the spontaneous emission spectrum in a four- 
level has been investigated by Ghafoor et al. 1 7J - This 
scheme was utilized further for localization of an atom 
during its motion on the classical standing- wave field . 
Thus observing a spontaneously emitted photon can lead 
to atom localization in a variety of systems. 

It is, however, important to note that from an experi- 
mental point of view, observation of a spontaneous emis- 
sion spectrum is very tricky and difficult. In this context, 
another scheme based on a three-level A-type system in- 
teracting with two fields — a probe laser field and a classi- 
cal standing- wave coupling field — has been used for atom 
localization by Paspalakis and Knight 19]. They observe 
that in the case of a weak probe field, measurement of 
the population in the upper level leads to subwavelength 
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localization of the atom during its motion in the standing 
wave. Thus in essence this scheme uses absorption of a 
probe field for atom localization. 

In this article, we describe another method for the 
localization of an atom in a standing-wave field based 
on electromagnetically induced transparency pcf. We 
consider a four-level system driven by two driving fields 
and a classical standing-wave field. A similar scheme 
was recently considered by us [2lJ to demonstrate phase 
and amplitude control of the group velocity of a weak 
probe field. The scheme is very similar to the one used 
by Ghafoor et al. for phase and amplitude control of 
the spontaneous emission 0] and atom localization . 
However, their consideration was based on monitoring 
the spontaneous emission properties of this scheme. In 
the present study we consider the absorption properties 
of the weak probe field by determining the susceptibility 
of the system at the probe frequency. We show that the 
probe absorption observed at appropriate frequencies lo- 
calizes the atom in the classical standing-wave field. We 
also investigate the effects of the amplitude and phase of 
the driving fields on the precision of localizing an atom. 
A novel feature observed by us as compared to the ab- 
sorption based scheme of Paspalakis and Knight is 
as follows: They observe four localization peaks in one 
wavelength of the standing- wave field — i.e., two each in 
the two half-wavelength regions. However, we obtain a 
more precise localization in the sense that for a certain 
choice of parameters we can confine the atom to one of 
the half-wavelength regions. We obtain two localization 
peaks that occur in one of the half-wavelength regions. 
We term this new domain of subwavelength localization 
as sub-half-wavelength localization. We would like to 
point out that this sub-half-wavelength localization has 
already been proposed by Ghafoor et al. , in the con- 
text of monitoring spontaneous emission spectrum. How- 
ever, they have not used the term sub-half-wavelength 
localization. 

This article is organized as follows. In Sec. [H] we intro- 
duce our model and give the basic equations and their 
solution to determine the susceptibility. In Sec. IIIII we 
study the behavior of the susceptibility along the nor- 
malized position coordinate along the standing wave for 
a variety of system parameters. This gives us conditions 
on the controllable parameters to attain localization of 
the atom as it passes through the standing-wave optical 
field in the cavity. We give an analytical expression for 
the appropriate detuning parameters to obtain sub-half- 
wavelength localization. Finally we present our conclu- 
sions in Sec. IIVI 



The internal energy level structure of the atom is shown 
in Fig.^b). The radiative decay rates from the levels \ai) 



II. MODEL AND EQUATIONS 

The schematics of the proposed scheme are shown in 
Fig. ^ We consider an atom, moving in the z direction, 
as it passes through a classical standing-wave field of a 
cavity. The cavity is taken to be aligned along the x axis. 
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FIG. 1: The Model: (a) The cavity supports the standing 
wave field (1) corresponding to Rabi frequency Qi . Two other 
fields (2, 3) are applied at an angle as shown. The atom enters 
the cavity along the z axis and interacts with the three drive 
fields. The whole process takes place in the x-z plane, (b) 
The energy level structure of the atom. Probe field, denoted 
by £ p , is detuned by an amount A from the |ai)-|c) transition. 
The fields (2, 3) shown in the (a) part of the figure correspond 
to the fields with Rabi frequencies Q2 and SI3, respectively. 
The decay rates from the upper levels |ai) and \a2) are taken 
to be 71 and 72, respectively. 

and \a,2) to level |c) are taken to be 71 and 72. The up- 
per level \a\) is coupled to the level |a2), and further the 
level 1 02} is coupled to level |6) via classical fields with 
Rabi frequencies and Sl 2 , respectively. In addition, 
the upper level \a\) is coupled to level \b) via a classical 
standing-wave field of frequency v and phase if, having 
Rabi frequency Oi. It should be noted that the Rabi fre- 
quency of the standing wave is position dependent and 
is taken to be £li(x) — f2isinra. Here £li(x) is defined 
to include the position dependence and K is the wave 
vector of the standing- wave field, defined as n = 2tt/X, 
where A is the wavelength of the standing-wave field of 
the cavity. We assume that the atom is initially in the 
state |c) and interacts with a weak probe field that is 
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near resonant with the |c) — > \a±) transition. The detun- 
ing of the probe on this transition is taken to be A. We 
assume that the center-of-mass position of the atom is 
nearly constant along the direction of the standing wave. 
Therefore we apply the Raman-Nath approximation and 
neglect the kinetic part of the atom from the Hamilto- 
nian p3 |. Under these circumstances, the Hamiltonian 
of the system in the rotating-wave approximation can be 
written as 



(1) 



where 



Ho = huj ai K) (ai\+hcj a2 \a 2 ) (a 2 |+?ia; b \b) {b\+Two c |c) (c| 

(2) 

and 



we only need following equations: 

Pa lC = -[*K 1C - Vp) + \li]~Pa lC + fae-^e tkxcos93 p a2C 
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Here are the frequencies of the states \i) and are the 
frequencies of the optical fields, and 8 2 and #3 are the an- 
gles made by the propagation direction of the fields fl 2 
and f2 3 with respect to the x axis, respectively. The sub- 
script p stands for the quantities corresponding to the 
probe field; i.e., £ p and v p are the amplitude and fre- 
quency of the probe field. Also p ai c is the dipole matrix 
element of the c) — > |ai) transition. For simplicity we 
assume that the Rabi frequencies f^i and Q 2 are real and 
O3 is complex — i.e., f2 3 = f2 3 e _I(p . This choice of impart- 
ing a carrying phase to field 3 is only for the convenience 
of calculations. As will become clear later, only the rel- 
ative phase of the three fields is important and absolute 
phases do not matter. The dynamics of the system is 
described using density matrix approach as 
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-{H,p}--{T,p} 



(4) 



where {T, p) — Yp + p]?. Here the decay rate is incorpo- 
rated into the equation by a relaxation matrix T, which is 
defined by the equation (n\T\m) — 7 rl i5 nTO . The detailed 
calculations of these equations are given in the Appendix. 

Our goal is to obtain information about the atomic po- 
sition from the susceptibility of the system at the probe 
frequency. Therefore we need to determine the steady 
state value of the off-diagonal density matrix element 
p aic . After necessary algebraic calculations and moving 
to appropriate rotating frames, we obtain a set of den- 
sity matrix equations [see Eq. Q24[l] . To determine p aiC 
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As we know, the dispersion and absorption are related 
to the susceptibility of the system and are determined by 
Paic- We take the probe field to be weak and calculate 
the polarization of the system to lowest order in £ p . We 
keep all the terms of the driving fields but keep only 
linear terms in the probe field. The atom is initially in 
the ground state |c); therefore, we use 



Pec — L > Pba x ~ U ' Pa 2 a t — U ' Paxax ~ u - \ D I 

Equations (JSJ can then be simplified considerably to ob- 
tain 

P aic = -(iA+ i 7 i)Pa lc + 5O3 e^e lkxcose *p a2C 

+ -ih sm kx p bc + i 1 
Pa 2C = -(iA+ l - l2 )p a2 c + 5O3 e^e-* kxcos9 *p ai c 

+ ^n 2 e tkxcos9 *p bC7 

1 

Pbc = ~i A pbc + ^ ^1 Sm KX Paic 



~(0) 
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Pa 2 c- 
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Here we have introduced the detuning of the probe field 
and the frequency difference between levels |ai) and |c), 

A = UJ ai c ~v v = u a2C + v 3 - V p = LU bc + V\ - v p . (8) 

We have also assumed 7 0C = 0. 

This set of equations can be solved analytically. Fol- 
lowing the treatment discussed in the Appendix, the 
off-diagonal density-matrix element corresponding to the 
probe transition is obtained as 



Pd\C Pd\C& 



1 



' - (fii - 4A 2 + 2t 72 A)£ pPaiC e-^ t , 

(9) 

where we have assumed 6* 3 = n/4, 6 2 = it/2 + it J A, and 
Y is given by 



Y = A + iB, 



(10) 
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with 

A = -8A 3 + 2A(ft 2 S i n 2 KX + q2 + ^ 

+ 2 7 i 72 A + Qin 2 03(e^ + e^ v ) sinra, 
B = 4A 2 ( 7l + 72 ) - 7l ft 2 - l2 Vi\ sin 2 kx. (11) 

The susceptibility at the probe frequency can be writ- 
ten as 



= 2Npa 1 cPaic e i Up t 



2N\p aiC \ 2 (ft 2 - 4A 2 + 2i 72 A) 



eo 



(12) 

where JV is the atom number density in the medium. The 
real and imaginary parts of susceptibility are given as 



X 



/ _ 2N\p aiC 



X 



eohZ 

u 2iV|p QlC | 2 



{(ft 2 -4A 2 )A + 2 72 AB}, (13) 
{2 72 AA-(ft 2 -4A 2 )£?)}, (14) 



where Z = YY* and x — x' + ix" ■ It is imperative 
to point out that the phase enters the susceptibility ex- 
pression only through the quantities A and Y. Even 
the phase dependence of Y is only through A. More- 
over, we observe that the phase-dependent term in A is 
ftift 2 ft3(e J(p + e~ lv ) sin kx. Thus the phase factor could 
very well have come from either of the three driving fields. 
As pointed out earlier, if all the fields had phase depen- 
dence, only the collective phase would be important and 



no individual phase-dependent terms would occur. This 
is because the Rabi frequencies ft^ in all the other terms 
appear through ft 2 , which is |fti| 2 for a complex Rabi 
frequency fij = |ftj|e l< ". The collective phase can be 
easily determined to be <f> — fa + fa — fa by repeating 
the susceptibility calculation. Here fa is the phase of the 
complex Rabi frequency ft^ of the zth driving field. 

In the next section we consider the imaginary part of 
the susceptibility x" m detail and obtain various condi- 
tions for subwavelength localization of the atom. 



III. RESULTS AND DISCUSSION 

We study expression (|14fl for the imaginary part of the 
susceptibility on the probe transition in greater detail in 
the following discussion. It is clear that X"— i.e., the 
probe absorption — depends on the controllable parame- 
ters of the system like probe field detuning and ampli- 
tudes and phases of the driving fields. 

Noting the dependence of x" on sinra, it is, in princi- 
ple, possible to obtain information about the x position 
of the atom as it passes through the cavity by measuring 
the probe absorption. Nevertheless, for precise localiza- 
tion of the atom the susceptibility should show maxima 
or peaks along the x coordinate. We obtain the condi- 
tions for the presence of peaks in x" m the discussion 
to follow. In the case of 72 = — i.e., the level |a 2 ) is 
metastable — Eq. (|14|) can be simplified as follows: 



„ _ 2N\p a J 2 7l (» 2 -4A 2 ) 2 

X he 1 2 (n 2 ~AA 2 ) 2 + [8A 3 -2A(fllsm 2 Kx + nl + nl)~2n 1 n 2 ^3COSipsmKx} 2 

2N\p aiC \ 2 7l (» 2 -4A 2 ) 2 

he 7l 2 (r> 2 ~4A 2 ) 2 +4A 2 r!f (sin K x-i?i) 2 (sinKa;- J R 2 ) 2 ' 1 ' 

where 

R Xft = — I -^2^3 cos ^± y / 0|Oi cos 2 </?-4A 2 [(ft 2 + ft 2 ) -4A 2 ]j . (16) 



From Eq. I|15[l it is clear that peaks will occur in x" a t 
x positions satisfying sin/tx = i?i. 2 . In other words, the 
probe absorption peaks at the spatial position defined by 



kx = sin 



1 



2Afti 



— £l 2 f23 cos ip 



± Jn 2 n 2 cos 2 ip - 4A 2 [(ft 2 + n 2 ) - 4A 2 ] 



± 717T, 



(17) 



where n is an integer. This leads to localization of 
atoms conditioned on probe absorption at a particular 



frequency corresponding to the value of A. 

We show the dependence of x" or probe absorption in 
arbitrary units versus the dimensionless x— coordinate in 
Fig. We show how the number of peaks, their posi- 
tions, and widths vary as the probe-field detuning, and 
driving field Rabi frequencies and the relative phase car- 
ried by the standing- wave field are changed. As seen from 
(a)-(c) in Fig. |21 the number of peaks is dependent on the 
detuning. As the detuning is decreased the peaks sepa- 
rate and we obtain four peaks in (b) and (c) compared 
to only two in (a). In (d) we show that, for a larger 
strength of the standing-wave cavity field, the peaks be- 
come sharp, leading to localization of the atom at one 
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FIG. 2: Dependence of the probe absorption on drive-field Rabi frequencies and detuning: plot of the imaginary part of the 
susceptibility in arbitrary units vs the dimensionless x coordinate kx along the standing wave in the cavity, kx runs from the 
values — 7r on the extreme left to n to the extreme right in each box. A vertical line is drawn at kx = to guide the eye. The 
common parameters are 0.2 = £la = fl — ji, 72 = 0, Hi = 371, and tp — n/2 unless specified otherwise, (a) A = 571, (6) 
A = I.471, (c) A = I.371, and (d) f2i = 2O71, A = 571. Notice that with the increase in the strength of the standing-wave field 
the peaks get sharper. Although the detuning is the same for (a) and (6), the number of peaks is different as the cavity-field 
strength is different. 



of the four possible positions. It can be noted that such 
localization is conditioned on the measurement of the ab- 
sorption of the probe field at a frequency corresponding 
to the chosen value of the detuning. 

The positions of the probe absorption maxima are 
strongly dependent on the probe-field frequency through 
its detuning A defined in Eq. JSJ. To clarify this point 



further we consider the relation sin/ix = R\,2 for various 
values of the phase parameter tp and solve for the detun- 
ing. This gives the values of the detunings as a function of 
kx to obtain probe absorption peaks, provided all other 
parameters are fixed. Taking O2 = = Q for simplicity 
we obtain 



A =4 



Oi sin kx zt y Q\ sin kx + 80 



°1,2 



A = ±~ ^/O 2 sin 2 kx + 20 2 = ± 5 { * /2) for tp = tt/2, 



or A = — — Qi sinra = 5^ for tp = 0, 



A = - [-Ox sin kx± y O 2 sin 2 kx + 8O 2 ] = 5$ 



or A = — 0\ sin kx = <5 3 for tp = tt . 



From these relations it can be seen that there is a strong 
phase dependence on the detuning values required to ob- 
tain probe absorption maxima. Also for each phase value 
there are several detuning values that would give rise to 
probe absorption maxima at each position. Namely, for 
tp = and 7T there are three different values given by 5^1 3 
and S 1^2 31 respectively. Similarly for tp = tt/2 detuning 
can take either of the values ±5^^ 2 \ 

To study this dependence in depth we plot the right- 



(18) 

(19) 
(20) 

I 

hand side (RHS) of the Eqs. (0-1(2011 at different posi- 
tions in one wavelength of the cavity standing- wave field 
in Fig. 13 These curves give the detuning value required 
to obtain peaks in the probe absorption at the corre- 
sponding position along the cavity axis. However, in 
an experiment a particular value for the detuning of the 
probe needs to be chosen to begin with. Therefore, for 
a given fixed detuning, represented by horizontal lines 
in the plots, the maximum conditions would be satis- 
fied only if the horizontal lines intersect with the curves 
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obtained from the solution of Eqs. lfT%|l - (|20)) . That is, 
there will be certain x positions at which probe absorp- 
tion maxima will occur. The curves obtained for ip = 
and ip = 7r happen to be related to each other by the 
transformation x — > —x. Moreover, there exists another 
symmetry for the phase values ip — 0, tt. The curves ob- 
tained by changing the sign of the detuning are the same 
as the original ones transformed according to x — > —x. 
Whereas, the curves for p = tt/2 are entirely different and 
the structure is independent of the sign of the detuning. 
Keeping these symmetries in mind we only consider pos- 
itive values of the detunings, corresponding to a probe 
frequency less than the probe transition frequency, for 
the rest of the discussion. The results for negative values 
of the detunings can be similarly obtained. 

It can be seen from the plots in Fig. [S] that the num- 
ber of intersecting points and their positions depend on 
the value of the detuning chosen as well as the relative 
phase ip of the drive fields. The points of intersection cor- 
respond to the probe absorption maxima. In Fig. [3] wc 
have considered several values of the detunings for differ- 
ent phase parameters. From the figure the positions and 
number of peaks to be observed in the probe field ab- 
sorption can be predicted. To understand this in detail 
we consider these special detuning values -labeled (a)-(l) 
in Fig. [3] and plot numerically evaluated x" m Fig- H 

It is straightforward to make contact with the results 
presented in Figs.^and^J For each value of the detuning, 
represented by horizontal lines in Fig. |3 we have a box 
showing the corresponding probe absorption curve along 
the cavity standing-wave field. The predictions about 
the number of peaks and their position for each detuning 
and phase value made by observing Fig.|2|can be verified 
by looking at the corresponding box in Fig. 01 Thus, 
we observe that the number of peaks and their positions 
depend on the phase and detuning values as long as other 
parameters are kept fixed. Using Fig. 0] we can make 
several observations. 

We obtain an interesting regime of atom localization 
for <p = and tt; the probe absorption peaks are situated 
in either nx — {— 7r, 0} or kx — {0,7r} half- wavelength 
regions along the cavity field, provided the detunings are 
chosen properly [see Fig. 0] boxes (b), (c), (d), (j), (k), 
and (I)]. The more interesting regime corresponds to the 
existence of just two peaks as opposed to the usual four. 
We have coined the term sub-half-wavelength localiza- 
tion for such an interesting regime of localization. It can 
be seen that for tp = tt/2, although the number of peaks 
and their widths vary as per the detuning value, there 
is no sub-half-wavelength localization possible. Neverth- 
less, the usual atom localization regime [See Fig. box 
(g)] is still available. Another interesting observation for 
the parameters ip = tt/2, A = [See (e)] gives no peaks in 
the probe absorption. The probe absorption is uniform 
over the length of the cavity, although the cavity field 
strength is different at different points. As seen from 
Fig. |3 this value of the detuning does not intersect with 
the solutions of the maximum condition sin^a; = Rx^. 



For all the results represented in this figure we have cho- 
sen positive values for the probe detuning; however, the 
results can be similarly obtained for negative detunings 
as well and they show similar characteristics. 

It is imperative to clarify the results obtained here. It 
can be seen that observing the peak in the absorption of 
the probe field of a particular frequency leads to local- 
ization of the atom. It means only if the atom passed 
through that position along the cavity field, one would 
obtain peak absorption for the probe. Another atom 
passing through the cavity may not pass through the po- 
sition corresponding to the absorption peak; in such a 
case, the probe field will not experience any absorption. 
This can be alternatively understood as the center-of- 
mass distribution of the atomic beam is not changed by 
the procedure considered here. It is only modulated by 
the coherent processes occurring in the system. There- 
fore, the localization of the atom is conditioned upon 
the observance of a peak in the absorption of the probe 
at particular frequency. However, the fact that a par- 
ticular frequency is associated with a particular location 
along the standing wave, one can envision the probe field 
carrying multiple frequencies. Such a configuration can 
then lead to localization of atoms at different locations 
conditioned upon the peak absorption of the correspond- 
ing frequency component of the probe. Such a scheme 
can be applied for atom lithography for generating ar- 
bitrary one-dimensional patterns. One can extend these 
arguments for generating arbitrary two dimensional pat- 
terns as well. The advantage our model provides is that 
the resolution is much larger compared to conventional 
subwavelength localization. We have succeeded in local- 
ization in the sub-half-wavelength domain. 

Thus, we have shown how to obtain sub- half- 
wavelength localization through monitoring the probe 
absorption at a particular frequency. We note that the 
atom is to be prepared in its ground state to start with 
[See Eq. JBJ] as opposed to schemes based on the ob- 
servation of the spontaneous emission spectrum (e.g., 
Ref. 0|) j where the atom needs to be prepared in its 
excited state. Thus the preparation stage is considerably 
simplified in our model. Moreover, as we need monitoring 
the probe absorption as opposed to spontaneous emission 
as in Ref. [l8j, we have a distinct advantage to offer as 
the absorption measurements are straightforward to real- 
ize in an experiment compared to the measurement of the 
spontaneous emission spectrum. In the following section 
we summarize our conclusions. 



IV. CONCLUSION 

We have presented a scheme for subwavelength lo- 
calization of a moving atom as it passes through the 
standing-wave field of a cavity. This allows us to de- 
termine the position of the moving atoms with high pre- 
cision as we show the possibility of confinement of the 
atom to a sub-half-wavelength region. Our scheme is 
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FIG. 3: Plot of the RHS of Eqs. 1181 - 12UH (in units of 71) vs the dimensionless x coordinate kx. Each curve gives the values 
of the detunings at given positions corresponding to peaks in the probe absorption. The common parameters are taken to be 
^2 = ^3 = = 2O71, 72 = 0, and fii = 3O71. Note that the character of these plots would change for different drive-field 
parameters. Interesting values of the detunings chosen for further analysis are shown by horizontal lines in each graph. The 
number of places at which the horizontal line intersects with the plotted curves corresponds to the number of peaks observed 
in the probe absorption, and the corresponding x coordinate gives the positions of the peaks in the cavity standing-wave field. 
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FIG. 4: Phase dependence of the localization for several values of the detuning parameter. Plot of the imaginary part of 
the susceptibility in arbitrary units vs the dimensionless x coordinate kx. kx runs from the values — ir on the extreme left to 
ir to the extreme right in each box. A vertical line is drawn at kx = to guide the eye. The common parameters are the 
same as in Fig.|3 The labels (a)-(l) correspond to the detuning values similarly labeled in Fig.|H] Notice sub-half-wavelength 
localization for boxes (b), (c), and (d) for phase ip = 0. Similarly, for ip = ir sub-half- wavelength localization exists as shown in 
boxes (j), (k), and (I). The term sub-half- wavelength localization is coined for a special regime of atom localization where the 
localization peaks are confined to either the range kx — {— ir,0} or kx — {0, ir}. It can be noted that for ip = ir/2 there is no 
sub-half-wavelength localization, although the number of peaks varies according to the detuning value chosen. The interesting 
results correspond to boxes (d) and (I) where there are only two peaks that are confined to a half-wavelength region on the 
cavity field. In box (e) the probe absorption is the same at all positions in the cavity field represented by the dark line appearing 
at the top boundary of the box. 



based on the measurement of absorption of a weak probe 
field by the atom. Measurement of the absorption of 
a weak probe field of prechosen frequency localizes the 
atom in real time. We have shown that the precision 
of the position measurement of the atom depends upon 



the amplitude and phase of the classical standing-wave 
field. The amplitude of standing-wave driving field when 
increased leads to line narrowing in the probe absorp- 
tion, thus giving increased precision in the position mea- 
surement whereas the phase of the standing- wave driving 
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field has important role in reducing the number of local- 
ization peaks from the usual four to two, leading to a new 
localization scheme which we call sub-half- wavelength lo- 
calization. Moreover, we show that the proper choice of 
probe frequency is very important in obtaining sub-half- 
wavelength localization. As our method is based on the 
measurement of the probe absorption, it has two dis- 
tinct advantages compared to similar methods based on 
observation of the spontaneous emission spectrum. Ab- 
sorption measurements are much easier to perform in a 
laboratory compared to monitoring of the spontaneous 
emission spectrum. Moreover, we do not require the 
atoms to be prepared in their excited states. We require 
that they are prepared in their ground state, which is a 
fairly routine task in atomic physics experiments. Thus 
the preparation stage is fairly straightforward. These ad- 
vantages suggest an easy experimental implementation of 
our scheme. 
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DETAILS OF THE SUSCEPTIBILITY 
CALCULATIONS 

Here we give details of the derivation of the density 
matrix equation and their solution for the calculation of 
the susceptibility. The density matrix equations are given 
by Eq.(4), so the (i, j)th element of the density matrix 
satisfies the equation 



Pij 



t {H ik , pkj - H k j ) - - y (Tik + Pik^kj ) ■ 

(21) 

Here the indices i,j, and k run over ai, 02, b and c . 
The decay matrix has only two nonzero elements given 
by r aic = 7! and T a2C = 7 2 . By utilizing the Hamilto- 
nian given in Eqs. (l)-(3) and Eq. 1)21(1. the off-diagonal 
density matrix elements can be shown to satisfy the fol- 
lowing set of equations: 

Pa ia2 = -(lUaias + 7;(7l +72))p Ql a 2 

_ io„ p -i<P p ikx cos 63 ivztt _ \ 

2 3 e e \Faiax Pa 2 a 2 ) 



Patb 



- io„ p -ikx cos 8 2 iv 2 t . i£pPaic -w p t 

2 2 Haib \ Pca 2 i 

= -(iu ai b + 7j7i)ah& 

- -tlx sinnxe~ Wlt (p aiai - pbb) 

1 (-) ikx cos 82 „—iv2t n 
r) i '2 c e Pa\a 2 



Pa 2 b = -(i^a 2 b + 7^2) Pa 



2% 



! n ikx cos 9 2 — iv 2 t l \ 
2^2 e e (Pa 2 a 2 - Pbb) 

+ ln 3 e i v e - ik3:cose °e il ' st p aib , 
Pai c = -(*w 0lC + \li)p aic + ^ 3 e-^e fccose3 e^Va 2C 
+ — fli sin/tx e~ Wlt pbc 

_ iSppa+c -iv p tt _ \ 

2^ e VPaiai Pec), 

Pa 2 c = - (We + \l2) P a 2 c + ^3 e~ ikx C ° S <V 3 < PaiC 
i 4 „ikxcos8 2 — iv 2 t 

-r 2 2 e e Pbc 

_ i£pPa lC -j Vp t 

2h Pa2ai ' 

Pbc = - (i u hc + 7 bc ) p bc + sin kx e lVlt p aic 

+ *_ Q2e -ikXCOS92 e iV2t pa2c _ !^l£ e -^t pbai . 

(22) 

The diagonal density matrix elements — i.e., the popu- 
lations of the atomic energy levels — can be determined in 
a similar manner. Here uiik corresponds to the energy dif- 
ference between the levels i, k, and Vk is the frequency of 
the fcth optical field. Next we transform these equations 
to appropriate rotating frame defined through 



Pca 2 

p ai b 

Pa 2 b 
Pa\a 2 
Pbc 
Pa\c 



i(v p -v 3 )t~ 
e Pca 2 , 

-i{V2+Vz)t ~ 

c f J a 1 a 2 , 

e- iv ^p aic . 



(23) 



The equations of motion for the density matrix ele- 
ments in the rotated frame take the following form: 



Pa\a 2 



i (^qiq 2 - v$) + ^(71 +72) 



Pan 



-fli sin kx e Wl p b , 



-f2i p -i<P „ikx cos 9 3 ( ~ _~ \ 

,"3 C e VPaiai Pa 2 a 2 ) 
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Pa 2 b 



-Oisin kx p ba2 

* O „~ikx cos9 2 ~ , i i£pPaic - 

2 S '2e p ai b H ^ — ^ ca 2' 

n 1 1- 

l (UJ ai b + ~7lJPai6 



N 1 1 - 



O ^zfcx COS 02 ( ~ ~ \ 



Pa%c 



i (u aiC - u p ) + -7i 



Pa%c 



Pa 2 c 



_0„ p-itfpikx cos 8 3 ~ 

-S2iSinKXp6c 2~ft 0°aiai ^Pcc), 

ik 2 c - (fp - ^3)] + ^72 P~a 2 c 



'Pa 2 a± i 



pbc 



2h 

i[lObc + [V\ - U p )] + 76c 



-n L sin kx p aic 



Pbc 



~Q 2 e 



-ikx cos 2 z i£pPaic - 



' Pa 2 c 



2h 



Pbat 



(24) 



Now we use the definition of the detuning from Eq. JSJ) 
into Eq. H24|) and after linearization with respect to the 
probe field we obtain Eq. 10 which is only a set of 
equations useful for determining the susceptibility of the 
medium at the probe frequency. 

Now we consider the details of the solution of Eq. J2J • 
This set of equations can be solved by writing them in 
the matrix form 



R{t) = —MR(t) + C , 



(25) 



where R(t) and C are the column vectors and M is a 
matrix: 



R — (Paic Pa 2 c Pbc) , 



C=(i^ 



M = 
( 



iA + - 7l ) 



.-Sl 3 e- i *e ikx 00803 --fiisinAtx 
2 2 



1 



--VL z e lv e- tkxcose3 (iA + - 72 



O „ikx cos 62 

--n 2 e 



V 



— ill sin kx 
2 



iA 



The formal solution of such equations is given by 

R(t) = f e-^-^Cdt' = M- X C. 

J — OO 

Finally by using Eq. (|27|1 . we obtain the solution, Eq. 



(26) 



(27) 
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